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Direct J/ij: and V' production rates at Tevatron are calculated in the kt- 
factorization approach within the color-singlet model. In this approach, the 
production rates are enhanced by a factor of 20 compared to the naive collinear 
parton model. However, the theoretical predictions are still below the exper- 
imental data by at least one order of magnitude. This means that to explain 
charmonium productions at Tevatron, we still need to call for the contribu- 
tions from color-octet channels or other production mechanisms. 
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Typeset using REVTgX 



Studies of heavy quarkonium production in high energy colhsions provide important 
information on both perturbative and nonperturbative QCD. In recent years, heavy quarko- 
nium production has attracted much attention from both theory and experiment. To explain 
the J/%1) and ip' surplus problem of large transverse momentum production at Tevatron [|l|, 
the color-octet production mechanism was introduced for the description of heavy quarko- 
nium production p|J^ based on the NRQCD factorization framework [Q. During the last 
few years, extensive studies have been performed for the test of this color-octet production 
mechanism 0. However, it was announced recently 0] that the color-octet contributions 
to charmonia productions at Tevatron are excluded by the experimental data under the 
/cj-factorization approach ||^,|[ , where the authors calculated Xcj production rates and found 
that the experimental data are successfully described only by the color-singlet contributions. 
Consequently, an urgent problem arises, which is whether the direct J/ip and ip' productions 
can also be described only by the color-singlet contributions in the kt factorization approach. 
In this paper, we will try to tackle this problem. We will calculate the color-singlet direct 
J/ip and ip' production in the kt factorization approach. We find the production rates are 
enhanced by a factor of about 20 compared to the conventional collinear parton model pre- 
dictions after considering the kt effects of the incident partons. However, the theoretical 
predictions are still below the experimental data by at least one order of magnitude. There- 
fore, to explain the Tevatron data on large transverse momentum charmonium productions, 
we still need to call for color-octet contributions or other production mechanisms. 

The /c(-factorization approach differs greatly from the conventional collinear approxima- 
tion because it takes the non-vanishing transverse momenta of the scattering partons into 
account. The conventional gluon densities are replaced by unintegrated gluon distributions 
which depend on the transverse momentum kt. 

The lowest order graphs for heavy quark pair production in kt factorization approach 
are plotted in Fig. 1, where we just give the typical graphs of The first two types of 
graphs (a) and (b) are conventional ones for the gluon-gluon fusion processes. The third 
type of graphs (c) are needed to preserve gauge invariance. To calculate heavy quarkonium 
production in this approach, we must project the QQ pair into a particular bound state. 
Because the QQ pairs in (b) and (c) are in color-octet, these two types of graphs will not 
contribute to heavy quarkonium production within the color-singlet model. In fact in 0, 
the productions of color-singlet P-wave charmonium states Xcj are calculated, where there 
are only type (a) graphs contributing to their productions. 

For the spin-triplet S'-wave charmonium states J /if} and ip' , in the lowest order, even 
type (a) graphs of Fig. 1 can not contribute to their productions because of C (charge) 
parity conservation. So, to produce the spin-triplet color-singlet S'-wave heavy quarkonium 
states, one must go to higher order (in terms of as) processes where an additional gluon 
must be emitted out. That is to say, in the lowest order, the spin-triplet S-wave quarkonia 
are only produced in the QQg production processes. In Fig. 2, we give the typical graphs 
for QQg production in the fc^-factorization approach, which include all contributing graphs 
in this order calculations in a gauge invariant formalism [Q. However, among these graphs 
only type (a) graphs of Fig. 2 contribute to S'-wave quarkonium productions, because other 
graphs violate C parity or color conservations. 

To calculate the production amplitude given by the type (a) graphs of Fig. 2, we make 
a Sudakov decomposition for every 4-momenta ki as 
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ki = aipi + (3iP2 + kiT, 



(1) 



where pi and p2 are the momenta of the incoming hadrons. In the high energy hmit, we 
have p\ = 0, P2 = 0; ^ind 2pi ■ p2 = s, where s is the cm. energy squared, and Pi are 
the momentum fractions of pi and p2 respectively. kiT is the transverse momentum, which 
satisfies 

kiT ■Pi = 0, kiT ■P2 = 0. (2) 
For the momenta of the incident gluons qi and q2, we have the following decomposition [§, 

qi = xipi + qiT, q2 = X2P2 + q2T- (3) 

That is to say, the longitudinal component of qi (^2) is only in the direction of light-like 
vector pi {P2). 

Using the above defined Sudakov variables, we can express the cross section for J/ip 
production from the type (a) graphs of Fig. 2 in the following form, 

da{pp J/ipX) = — — d pTd qiTd q2T 



64 X 16(27r)4 qfj, q^^ qiTq^T 

(4) 

where pt is the transverse momentum of J/ip., and are the momenta fractions of pi 
carried by J/ip and the outgoing gluon. The amplitude Aq describes J/?/' production in the 
gluon-gluon fusion processes g + g J/ip + g, with the off-shellness of the two incident 
gluons being q\ = —qf^ and g| = — respectively. We have checked that the amplitude 
Ao vanishes when qiT — > or g2T —>■ 0, which is required by the gauge invariance of this 
approach. 

/(x; q^) is the unintegrated gluon distribution, which is related to the conventional gluon 
distribution by 

2\ divq 



xg{x,fi)= -j^f{x-kT). (5) 

J rvrp 

The unintegrated gluon distribution /(x; q^) includes the evolution in x and k^ by the BFKL 
and DGLAP equations. In the nonperturbative region of small k^ the unintegrated gluon 



distribution is not known, so in practice the above equation may be rewritten as [H lO 



xg{x,f/) = xg{x,qo) + -r^f{x]kl), (6) 

which introduces a priori unknown initial scale q^ and the initial gluon distribution xg{x, q^). 
According to this procedure, the amplitude has the following decomposition form |p|. 



S{qiT,q2T) 



\Ao{qiT,q2T) 

9'it9'2T 



= 5(0, 0)e{ql - qlT)e{ql - qj^) + S{q,T, 0)^(g?^ - ql)e{ql - ql^) 

+ 5(0, q2T)e{ql - qlT)0{qlT - ql) + S{q,T, q2T)e{qlT - ql)0{qlT - ql)- (7) 
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The function S here is similar to the function J of P] called the impact factor, which is, 
in effect, an off-shell but gauge invariant cross section [§]. We note that 5'(0, 0) represents 
the cross section for gluon-gluon fusion processes with no kt effects associated with the two 
incident gluons. Substituting this part of S into Eq. @, we can reproduce the conventional 
collinear parton model results for S-wave quarkonium production. 



For numerical calculations, we choose the unintegrated gluon distribution of 11 , where 



the authors determined it by using a combination of DGLAP and BFKL equations. With 
the initial conditions 

ql = lGeV\ xgix, q^) = 1.57(1 - x^-^ 

they obtained an excellent fit to F2{x, Q^) data over a wide range of x and . As in Refs. 
||1^,|6[|, we set the scales iJ? for strong coupling constant Q;s(/i^) in the amplitude squared 
|y4op to be q\j. for the interaction vertex associated with the incident gluon gi, g^T the 
vertex associated with and p\ + (m = 2mc is the quarkonium mass) for the vertex 
associated with the outgoing gluon attached to the cc line. 

In the following we present the numerical results for color-singlet direct S-wave charmonia 
productions at Tevatron, compared with the experimental data [|l|. Fig. 3 is for J/ip, and 
Fig. 4 is for %p' . In the calculations, we use charm quark mass as ttlc = 1.5GeV, and the 
color-singlet matrix elements as (Of(^Si)) = l.OSGeV^ and {Of'i^Si)) = Q.67GeV^ 0. 
To compare with the experimental data we have imposed a cut on the rapidity of J/ip and 
ijj': |r]| < 0.6. From these two figures, we can see that the production rates are enhanced by 
a factor of about 20 compared to the conventional collinear parton model predictions after 
considering the kt effects of the incident gluons. However, the theoretical predictions are 
still below the experimental data by at least one order of magnitude. Therefore, to explain 
the Tevatron data on large transverse momentum charmonium productions, we still need to 
call for color-octet contributions or other production mechanisms. 

In conclusion, we have calculated direct J/ip and ijj' production at Tevatron in kt- 
factorization approach within the color-singlet model. We found that the production rates 
are enhanced by a factor of 20 compared to the conventional collinear parton model. How- 
ever, the theoretical predictions are still below the experimental data by at least one order 
of magnitude. This means that to explain charmonium productions at Tevatron, we still 
need to call for the contributions from color-octet channels or other production mechanisms. 
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Figure Captions 



FIG. 1. The typical graphs for QQ production in kt factorization approach in gauge invariant 
formahsm. 

FIG. 2. The typical graphs for QQg production in kt factorization approach in gauge 
invariant formalism. 

FIG. 3. The color-singlet predictions for large pt direct J/ if] production compared with the 
experimental data from Tevatron. The dashed hne denotes the result in the conventional 
coUinear parton model, and the solid line is the result in the kt factorization approach. 
FIG. 4. The same as Fig. 3 but for i^' . 
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FIGURES 
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Fig.l The typical graplis for QQ production in fej factorization approach in 
gauge invariant formalism. 
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